Theory of Optical Orientation in n-type Semiconductors 
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Time resolved measurements of magnetization in n-GaAs have revealed a rich array of spin de- 
coherence processes, and have shown that fairly long lifetimes (~ 100 ns) can be achieved under 
certain circumstances. In time-resolved Faraday rotation and time-resolved Kerr rotation the evo- 
lution of the magnetization can be followed as a function of temperature, applied field, doping level 
and excitation level. We present a theory for the spin relaxation in n-GaAs based on a set of rate 
equations for two interacting thermalized subsystems of spins: localized states on donor sites and 
itinerant states in the conduction band. The conduction band spins relax by scattering from defects 
or phonons through the D'yakonov-Perel' mechanism, while the localized spins relax by interacting 
with phonons (when in an applied field) or through the Dzyaloshinskii-Moriya interaction. In this 
model, numerous features of the data, including puzzling temperature and doping dependences of 
the relaxation time, find an explanation. 
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Spin coherence in semiconductors is attracting renewed 
attention due to the prospects of spintronics - informa- 
tion storage and processing using spin rather than charge 
degrees of freedom, and by the idea that spins in semi- 
conductors could serve as qubits for quantum comput- 
ers 0,|3|- For qubit appHcations, the spin degreees of 
freedom must be coherent, which necessitates a detailed 
understanding of the processes that limit spin lifetimes. 
Time-resolved measurements on n-type systems have re- 
vealed an array of decoherence processes, and have shown 
that fairly long lifetimes (greater than 100 ns) can be 
achieved in n-GaAs[3,0- In time-resolved Faraday rota- 
tion and time-resolved Kerr rotation the evolution of the 
magnetization can be followed as a function of temper- 
ature, applied field, doping level, and the intensity and 
duration of the pump pulse. Results on different materi- 
als (GaAs, GaN, ZnSe) are similar, pointing to universal- 
ity in the phenomena. Awschalom and Samarth 0| have 
reviewed the experimental situation. 

Our theory provides a systematic framework for in- 
vestigating the wide range of parameters studied in op- 
tical orientation experiments on n-type semiconductors. 
Previous theory has concentrated on either higher tem- 
peratures 3|, or on very low temperatures and very low 
magnetic fields Q- We account for certain puzzling ex- 
perimental observations by having two distinct types of 
spin states: localized donor states and itinerant conduc- 
tion band states with characteristic spin relaxation rates 
1/t; and 1/tc, respectively. In addition, a fast cross- 
relaxation rate, 1/tc,., between the localized and itiner- 
ant spins is a crucial feature of these systems, leading 
to the largest relaxation rate, either 1/tc or 1/r/, usually 
tending to dominate the spin dynamics for the whole sys- 
tem. Examples of this behavior are shown in Figs. 1 and 
2 below. In this paper we focus on presenting the the- 
ory and applying it to n-GaAs, reserving more extensive 



comparision to experiment for a later publication. 

The cross-relaxation between localized and itinerant 
spins occurs by the usual exchange interaction Hamilto- 



J 



i,k 



where the sum runs over impurity spins i and conduc- 
tion band states k and V is the volume of the system. 
This Hamiltonian conserves total spin and cannot, by it- 
self, relax the magnetization. However, it can transfer 
spin from localized to itinerant states. J may be esti- 
mated as J ~ ~e^a%. Here qb = 10.4 nm is the effective 
Bohr radius for an impurity. If virtual excitations to the 
upper Hubbard band are important, this estimate could 
be reduced and J could even change sign. The sign is 
actually not important for our purposes, since all exper- 
imental temperatures are well above the Kondo temper- 
ature. For n-GaAs an order-of-magnitude estimate for 
the cross-relaxation rate is l/rcr ^ (1 ps~^)(nimp/no), 
where uq — 10^^ cm~^ is a fiducial density. 

The spin orientation is created by a circularly- 
polarized optical pump pulse about 100 fs long tuned 
near the band gap, creating particle-hole pairs. The va- 
lence band holes depolarize quickly and fast recombina- 
tion (on a time scale of 50-100 ps) leaves the conduction 
band and localized donor state system with a net spin 
polarization along the propagation direction of the beam 
(z-direction). The time evolution of this polarization 
is tracked by applying a transverse magnetic field in the 
x-direction (Voigt geometry). The resulting precession 
about the x-axis and concomitant decay are measured op- 
tically, with 1 /T2 the relaxation rate of the macroscopic 
transverse magnetization. 

Our theory may be formalized by writing modified 
Bloch equations for the magnetization for times after 
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recombination (t > 100 ps after the end of the pump 
pulse). The holes have recombined and spin-conserving 
processes have thermalized the system subject to the 
constraint that the magnetization retains the polariza- 
tion produced by the initial excitation process. There 
are then two thermalized subsystems of electrons at am- 
bient temperature with relative occupations determined 
by standard thermodynamic methods. The localization 
of conduction band electrons onto impurity sites takes 
place at a temperature scale Timp determined by riimp- 
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50 K. Denote 



the localized and conduction band densities by ni and 
ric, with riimp = ni + ric- In the experiments, N^x, the 
density of electrons excited by the pump pulse, is small, 
Nex ^ nimp, except possibly for nominally insulating 
samples, which we discuss briefly below. 

We work in the frame which rotates about the a;-axis 
at a rate g* fisB /h, and is along the z-axis at time t — 0. 
In this frame, the dynamics are governed by 
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where -I- and — denote up and down spins in the rotating 
frame and 7cr = n^Tcr- By rewriting Eq. 1-4 in terms 
of the total densities, ni = ni+ -\- ni^ and ric = nc+ + 
ric- and the magnetization densities, m/ — ni^ — ni^ 



and rric 
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we find that the total densities are 



time independent, dni/dt = and dric/dt 
magnetization densities are determined by 
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and the 
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In general, the time dependence of the total magneti- 
zation m{t) = rndt) + mi(t) is a sum of two exponentials, 
exp (— r+t) and exp (—T-t) (behavior observed in exper- 
iments on n-GaN|6j) with eigenvalues 
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where S is given by 
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For n-GaAs we are in the regime I/tc- ^ l/^c 
1/t;, where the eigenvalues give two very different re- 
laxation rates: a very rapid relaxation given by r_|- k, 
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FIG. 1: Plot of IjT^ versus temperature. The data are 
from Ref. 0, with solid dots for B = T and solid squares 
for _B = 4 T. The lines connecting the data points are 
guides for the eye. The curves are a least-squares fit of 
Eq. 5 to the data. Dashed-dotted curve: (ni/nimp)(l/rDAf), 
Dashed curve: (n;/nimp)(l/T"s-ph) for _B = 4 T, Dotted curve: 
(nc/ni,„p)(l/rc), Solid curves: total XjT^. For B = T, 



1/r. 



-ph 



0. Inset: Tp versus temperature for a: 



10"" cm-' and b: 



o/{noTcr), with a timescale on the order of picosec- 



onds and a slower relaxation given by 
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with a timescale on the order of tens of nanoseconds. 
Given expressions for 1/t; and 1/tc, Eq. 9 gives the cal- 
culated total relaxation rate from our theory. This is 
the appropriate quantity to compare to the single expo- 
nential time dependence observed in experiments 
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There are various processes that can relax the non- 
equilibrium magnetization produced in optical orienta- 
tion experiments. The conduction band processes have 
been well studied 7], while the relaxation mechanisms for 
localized electrons are less well understood. 

D'yakonov-Perel' Mechanism Conduction band 
electron spins in n-GaAs relax primarily by the 
D'yakonov-Perel' (DP) mechanism due to lack of in- 
version symmetry in III-V systems. Lack of inversion 
symmetry, together with spin-orbit coupling, gives an ef- 
fective ^-dependent magnetic field, causing the spin of an 
itinerant electron to precess about an axis related to k. 
The precession frequency for an electron at wavenum- 
ber k is flopik), and the DP relaxation comes from 
switching the precession axis by scattering from one k- 
vector to another. Assuming s-wave scattering, one ob- 
tains 1/T£)p(fc) = 2n'jjp{k)Tp{k)/3, where Tp{k) is the 
momentum relaxation time. Averaging this expression 
over the Boltzmann distribution, and using the results 
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FIG. 2: Plot of 1/To versus applied magnetic field. The 
data are from Ref. |3, with solid dots for Uimp = 10^^ cm~^ 
and solid squares for rUmp = 10^* cm~^. The lines connect- 
ing the data points are guides for the eye. The curves are 
a least-squares fit of Eq. 5 to the data. Dashed-dotted 
curves: {ni/nimp)(l/TDM) for the two densities, Dashed 
curve: {ni/ni,np{l/Ts^ph), Solid curves: total 1/T2. For 
T = 5 K, (nc/n,„j,)(l/r,) < 0.001 ns-\ 

of Fishman and Lampel for the momentum average 0, 
we find a spin relaxation rate 1/tjjp — aopT^Tp, with 
a_Dp(th) = 9.0 X 10^^" K^^ps^^. Here Tp is the aver- 
age momentum relaxation time which has a complicated 
temperature and doping dependence (shown in Fig. 1) 
best taken from mobility data, /ie = eTp/m* [l^. 

Elliot- Yafet Mechanism Conduction band electron 
spins can also relax via the EUiot- Yafet (EY) |T3 mech- 
anism, due to ordinary impurity scattering from state 
k to state k' . With spin-orbit coupling the initial and 
final eigenstates are not eigenstates of S'^, the spin pro- 
jection operator, so this process relaxes the spin. One 
finds \/tey = aeyT^/rp, where a_Ey(tli) = 8 x 10^^° 
K~^. For the experimental parameters discussed below 
{nc/nimp){^/TEY) is three orders of magnitude smaller 
than the leading contributions to l/Tl^ . Thus, to simplify 
our analysis of experiments on n-GaAs, we set asY = 
and 1/tc = 1/tdp- 

Spins localized on donor sites cannot relax by the 
same scattering dependent processes that relax conduc- 
tion band spins. As a rule, relaxation times for localized 
states are longer than for itinerant states, due to phase 
space effects. Lifetimes for localized states can be very 
long: times in excess of 10'^ s have been measured for 
donor bound states of phosphorus-doped silicon [l^ . 

Spin-phonon Mechanism Acoustic phonons can re- 
lax localized spins by dephasing, due to spin-orbit cou- 
pling mixing spin states, if an external field B breaks the 
time-reversal symmetry present at zero field (Van Vleck 
cancellation). With B in the z— direction, the relevant 
term in the Hamiltonian for a single spin takes the form 

Hs-ph = CfiBBa.^A,, (10) 



where is the dilatation at a donor site i and C is 
a constant. The dilatation modulates g* , the effec- 
tive factor, which is given by g* = 2[1 — (m* /m — 
l)Aso/{3Eg + 2Aso)], where Eg = 1.4 eV is the en- 
ergy gap, m*/m = 0.067 is the ratio of the effective 
mass to the bare mass, and Ago — 0.344 eV is the 
spin-orbit splitting of the valence bands 0. We also 
have that Eg{Ai) = Eg{0) - {9eV)Ai and the effective 
mass depends on the gap as to* ~ ^/Eg. These facts 
allow us to estimate « 550. There is no gener- 
ally accepted theory of the multi-spin relaxation rate 
I/T2* that results from such a Hamiltonian (the single- 
spin rate I/T2 has recently been calculated in silicon 
[iJI). We can obtain a simple estimate using Redfield 
theory, which yields l/rg-ph = C^JS^ (A^) Tp/i, where 
Tpii is the phonon correlation time 'l5'| and A is the av- 
erage dilatation for the occupied donors. This leads 
finally to l/r,_p,, = a,^phB^T^f{T), where f[T) = 

JH"^'^ [(1/2) + (e^ - l)-i] dx in a Debye model for 
the phonons and 6d = 343 K for GaAs. We note that 
this theory is not likely to be valid at higher tempera- 
tures, where multi-phonon and Orbach processes become 
important. This issue is not settled even in insulators, 
and we defer full consideration of it to a later publication. 

Nuclear Hyperfine Mechanism A localized elec- 
tron spin can relax through the hyperfine interaction with 
the N « 10^ nuclei with which it is in contact. The nu- 
clei are randomly oriented under most conditions, and 
the associated field felt by the electron is A/^/NgfiB ~ 
lQ-2 -vvhere A is the hyperfine constant. The corre- 
sponding precession frequency is ujn ~ 10^ s~^. There 
is also a characteristic time for the nuclei Tn2 ~ 10"'* s, 
which comes from the nuclear dipole-dipole interaction. 
Hence uJNTn2 « 10^, and we are in the regime where the 
effective random field fluctuates slowly compared to the 
precession of the spin. The relaxation time from coupling 
to the nuclei, 1/T„tjc, is not expected to have strong tem- 
perature dependence in the range T > 1 K considered 
here and there should be no field dependence as long as 
B > 10^^ T. We treat 1/Tnuc as a constant. 

Dzyaloshinskii-Moriya Mechanism Localized elec- 
tron spins can relax by the Dzyaloshinskii-Moriya (DM) 
interaction 16]. This interaction, arising from spin-orbit 
coupling, produces a term proportional to b-siXS2 where 
b is related to the interspin separation, and to the ex- 
change integral between the wavefunction on sites 1 and 
2. This interaction is not isotropic in spin space and 
can therefore relax the spins. The calculation of the 
effect of this term on l/Tj* is not straightforward, since 
it involves aspects of the spin glass problem that are not 
entirely solved. Gor'kov and Krotkov have given the 
first term in a density expansion. We use their expres- 
sion, though with a more general distance dependence 
for the exchange interaction as a first step toward 
a theory valid at higher impurity densities. We find 
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1/tda/ = aoMriimpa^fDMiT). Here aDM(th) = 0.01 
ns^^ and the weakly temperature dependent, dimension- 
less function foMiT) « 32 at T = 5 K. 

The total relaxation rate for the localized spins in our 
theory is given by l/n = l/Ts-ph + l/r„„c + ^/tdm- 

In Figs. 1 and 2 we compare the results of our theory 
to experimental data on n-GaAs at riimp = 10^^ cm~^ 
and riimp = 10^^ cm~'^. Our procedure is as follows. 
Each mechanism above has very definite field, temper- 
ature, and doping dependence. The overall constant 
factor for each is less certain. The data are consistent 
with 1/t„„c = 0. However, for reasons explained be- 
low, the data do not set tight limits on 1/t„„c- We 
do a least-squares fit to the complete data set shown 
in Figs. 1 and 2 using the three remaining adjustable 
parameters, with the results aDp(exp) = 8.6 x 10^^° 
K-3ps-2, as-ph(exp) = 2.2 x 10"" T-^R-^ns"!, and 
Q!_DM(exp) = 0.031 ns~^. These values are in satisfac- 
tory agreement with a^p (th) = 9.0x10"^" K~^ps~^ and 
a_DM(th) = 0.01 ns~^. The remaining value as-p/i(exp) 
is best viewed as an estimate of the phonon correlation 
time, Tph = 1.2 x 10~^ ns. This seems reasonable, given 
that the inverse Debye frequency Ti/kBdo = 2.2 x 10~^ 
ns. But true comparison of theory and experiment for 
this prefactor awaits a more comprehensive theory of the 
phonon relaxation, as noted above. 

The fits against temperature at riimp = 10^^ cm~'^ for 
_B = T and B = 4 T are shown in Fig. 1. There are two 
surprising points about the data: (1) I/T2* is independent 
of i? at high T; (2) the T dependence is non-monotonic at 
higher fields. Point (1) is explained by noting that in our 
theory all B dependence comes from I/ts-ph- For this 
doping, the localized states are completely depopulated 
at high T. Point (2) is more subtle. In our theory, 
though both 1/tc and 1/t; increase with increasing T, 
their contributions to l/Tj* are weighted by nd riimp and 
ni/ riimp, respectively. Starting at low T, where localized 
states dominate, we have I/T2 decreasing with increasing 
T due to the decrease of ni /riimp as T increases. After 
the localized states are depopulated 1/tc dominates and 
I/T2 is an increasing function of T. 

The fits against field for riimp = 10^® cm^"^ and 
riimp = 10^^ cm~'^ at T = 5 K are shown in Fig. 2. These 
results are also surprising: there is strong enhancement 
of I/T2 by B at low doping, while at high doping, the 
dependence is quite weak. In this case, the explanation 
relies on the DM contribution. At high doping 1/tdj\/ 
dominates because of the short impurity-impurity spac- 
ing and consequent fast relaxation. This rate is field- 
independent. At low doping l/Ts-pu with its strong B 
dependence is more important. 

Experiments 0,01 show riimp = 10^^ cm~'^ is the "op- 
timal" (smallest l/Tj^*) doping value at low T. Going to 
lower dopings (data not shown) increases l/Tj*. This is 
due to a combination of effects. Increasing the doping 
at first raises the number of localized states, which have 



intrinsically much smaller decay rates. However, this 
process stops when the impurity states begin to overlap, 
and 1/tdm dominates. This is why the minimum l/Tj 
is near the metal-insulator transition. These conclusions 
are consistent with those of Ref. Q at small B and T. 

From an examination of Fig. 2, we can see why adding 
I/tViuc does not significantly improve or worsen the fit. 
Doing so adds a constant to I/T2, which would move 
all theory curves rigidly upward, not much affecting the 
overall goodness of fit. Thus we cannot get a meaningful 
limit on 1/t„„c with this data set. The question of how 
to average over nuclear degrees of freedom to find 1 / r„„c 
is not completely clear, and there are different results in 
the recent literature 0, [20j . 

The main shortcoming in our theory lies in the high 
doping regime, where spin-glass effects become impor- 
tant. This is shown by the B = Q points in Fig. 2, where 
I/T2 is dominated by I/tdm- Since we use a low-density 
(pairwise correlations only) expansion for 1/tdm, the ra- 
tio of the theory points is exactly the ratio between the 
densities, which is too small compared with experiment. 
A better theory would approach freezing as a collective, 
not pairwise, effect. Apart from this, all the qualitative 
features of the data are explained in our picture. 

Our analysis shows that the complicated B and T de- 
pendences for I/T2 observed in experiments are due to 
having two strongly interacting subsystems of spins: one 
localized and the other itinerant. The coupling of local- 
ized spins to phonons gives rise to the unusual magnetic 
field depndences of the relaxation rates. Overall, the 
very different T and B dependences for I/tc and 1/r/ 
coupled with population effects give the wide range of 
experimental phenomena observed. 
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